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Abstract 

The method of derivation of equations describing the evolution of the neutral Bose particle system 
at low temperatures is described. Despite the fact that we consider the neutral particles we account 
the short-range interaction between particles. As important limiting case we consider the particles 
in the Bose-Einstein condensate (BEC) state. This method is called the method of quantum hy- 
drodynamics, because natural for of the quantum mechanics rewritten in terms of material fields 
of observable quantities in three dimensional space is the set of equations, which look like the hy- 
drodynamics equations. It can be shown that from the quantum hydrodynamics (QHD) equations 
can be derived macroscopic non-linear Schrodinger equation. Most famous non-linear Schrodinger 
equation is the Gross-Pitaevskii (GP) equation, which is non-linear Schrodinger equation with non- 
linearities of the third degree. Non-linear Schrodinger equation defines the wave function in the 
medium or order parameter, which is a macroscopic parameter. There are generalizations of the 
GP equation. New term appears in the GP equation at account of the three-particle interaction. 
This term contains nonlinearity of the fifth degree. At more detailed, in comparison with the GP 
equation, account of the two particle interaction we come to the non-local non-linear Schrodinger 
equation. This equation contains spatial derivatives of the order parameter in the non-linear terms 
caused by the interaction. Particularly, non-local non-linear Schrodinger equation arises at the 
consideration of the two-particle interaction up to the third order by the interaction radius. In this 
terminology the GP equation corresponds to the first order by the interaction radius. For the BEC 
of the neutral particles with anisotropic long-range dipole-dipole interaction the generalization of 
the GP equation was also suggested. Detailed analyses of the applicability conditions shows that 
this equation valid for the system of dipoles parallel to each other, which do not change their 
direction, and where the dipole-dipole interaction interferences motion of particles. All described 
non-linear Schrodinger equation can be derived from the corresponding, and in some cases more 
general, QHD equations. And, for all described cases and types of interaction, the equations of the 
QHD can be derived directly from the microscopic many-particle Schrodinger equation. This chap- 
ter is dedicated to the description of the method of the QHD. We show the method of derivation of 
the hydrodynamics equation and present the method of derivation of the corresponding non-linear 
Schrodinger equation from the QHD equations. During derivation we admit advantages of the 
QHD method. The QHD equations contain information about thermal motion of the particles, 
when we consider the BEC evolution we have to neglect by the contribution of the thermal motion. 



2 



In the chapter we illustrate contribution of the temperature on the Bose particle dynamics at the 
example of the system of Bose particles with the three particle interaction in the first order by the 
interaction radius approximation. Finally, we discuss new properties of the dispersion dependency 
of the eigenwaves in the BEC obtained recently by means of the QHD. 

PACS numbers: 
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I. INTRODUCTION 



The model of a slightly non-ideal Bose gas at almost zero temperature is a fundamental 
model for the BEC description giving Bose-type energy spectrum suggested by N.N. Bo- 
goliubov The non-linear Schrodinger equation for the wave function of BEC was 

suggested later, which is called the Gross-Pitaevskii equation . It is a great tool for 

inhomogeneous BEC studying. This equation contains nonlinearity of third degree. The 
Gross-Pitaevskii equation was generalized for the BEC of particles with three-particle inter- 
action {4J], (see also more recent papers (g]— js] ) . In this case a nonlinearity of fifth degree is 
added to describe three-particle interaction. A generalization of the Gross-Pitaevskii equa- 



tion for dipolar BEC was also suggested 



13] . It is well-known that the Gross-Pitaevskii 



equation might be presented in hydro dynamical form, analogously to the hydrodynamical 
representation of the one-particle Schrodinger equation. 

We consider evolution of the BEC from point of view of the quantum hydro dynamic 
method. Initially the QHD method was developed for quantum plasma jl^]-jl3]. Quantum 
hydrodynamics of ultracold quantum gases was obtained later [ji]]. The QHD method 
allows us to study evolution of spinning particles 15], 16) and particles having electric 
dipole moment 17]. Method of the QHD allows us to derive equations of the quantum 
observable directly from many-particle Schrodinger equation. These equations arise as a 
chain of equations, and the chain has to be truncated. A way of truncation depends on 
the kind of interparticle interaction and a set of physical process we want to consider. The 
approximation of the self-consistent field approximation gives us possibility to get a closed set 
of equation for particles with a long-range interaction (we use it for modeling of the dipole- 
dipole interaction). A short-range interaction requires another treatment. The procedure 
was developed in Ref . 18] , and it has no specific name. We will present detailed description 
of this method in this chapter below. 



II. GENERAL STRUCTURE OF THE QHD 

This chapter is dedicated to description of the quantum hydrodynamics method in appli- 
cation to the dynamics of Bose- Einstein condensate (BEC) of neutral atoms 18] . Certainly, 



we will describe relation between the quantum hydrodynamics and the Gross-Pitaevskii 
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equation, which is usually used for description of the BEC evolution, but we start from ba- 
sic ideas lay behind the quantum hydrodynamics, which actually follow from basic principles 
of quantum mechanics. One of the main goals we want to show is that the quantum hydro- 
dynamics is a natural way of description of quantum many-particle phenomenon, along with 
the second quantization. Then we study dynamics of a quantum physical quantity we have 
deal with the operator of this variable. Knowledge of wave function allows us to calculate 



quantum mechanical average of the physical quantity [19] 

<L>= ^*L^dR, 



then the wave function satisfy to the Schrodinger equation 

k = E 4^ §cA + ? Kxi(r " t] + ^ ? U{1 Ti ~ r ' l} ' (1) 

where pf = —1T1V* is the momentum operator of the i-th particle, m ; is the mass of the i-th 
particle, and = U(\ r$ — r 3 - |) is the interaction potential. 

If we study a many-particle system we should choose useful physical quantities and find 
operators for these quantities. Development of the QHD shows that the first of the useful 
variables is the particle concentration, and it is easy to find it's operator in the coordinate 



representation. We just need to take the microscopic concentration of particles [20[, 21] 

N 
i=l 

Using this operator we get the quantum mechanical definition of the particle concentration 

'>*) = / dRj2^-^ + (RM(R,t), (2) 



n r, 



where dR = YliLi ^ r «- 

We can use the Schrodinger equation to study evolution of the concentration. For this 
purpose we differentiate the definition of concentration and obtain the continuity equation 

d t n + Vj = 0, (3) 

where 

r(r,t) = J cLR^5(r-r^(^) + ^ (4) 



is the momentum density. 

Analogous derivation, but for one-particle Schrodinger equation, is presented in many 
books in quantum physics. In one particle case our definition (J25 gives well-known formula 
for density of probability tp*if). 

Our definition (J2J) has one more important feature. Many-particle wave function ip(R, t) is 
defined in 3N dimensional configuration space, whereas many-particle phenomenon realize in 
three dimensional physical space. Thus, if we want to have a description method presenting 
this evolution in natural way we should make projection of 3N dimensional evolution of 
psi function on three dimensional physical space. However, we do not need to solve this 
problem since we have solved this problem by using of the definition ([2]) for the particle 
concentration. 

Having definition of the particle concentration we can derive the set of the QHD equa- 
tion from the many-particle Schrodinger equation. We have already derived the continuity 
equation (j3J), this equation contains new physical quantity. It is the particles current or the 
momentum density. To continue derivation of the QHD set of equations we need to derive an 
evolution equation of the momentum density, which called the momentum balance equation 
or the Euler equation. For this purpose we need to differentiate the momentum density with 
respect to time and use the Schrodinger equation for time derivatives of the wave function. 
In the result we find 

d t f + d p Ii aP = F a , (5) 

where U al3 is the momentum flux tensor, and F a is the force density. 

We have found general form of the quantum Euler equation. The momentum flux tensor 
describes kinetic and kinematic properties. Kinematic properties of spinning and spinless 
particles are different. In this chapter we will consider spinless particles and we will consider 
general properties and explicit form of the momentum flux tensor. Whereas, the force 
density is determined by an explicit form of the interaction in the basic Hamiltonian. We 
will consider two- and three-particle short range interaction and the long-range dipole-dipole 
interaction of particles having electric dipole moment. For simplicity we start derivation of 
the QHD equation for particles with the two-particle short-range interaction. Nevertheless 
it is enough to derive the Gross-Pitaevskii equation. 

We have described general picture. Now we present derivation of the explicit form of the 
momentum flux tensor and it's connection with the velocity field. 
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Explicit form of the momentum flux tensor arising in the Euler equation (jSJ) is 
ir*(r,f) = J dR^div - ri )^-(^ + (R,t)(ptp^)(R,t) 

+ (p?^) + (R,t)(pfy)(R,t) + c.c^. (6) 

To present the velocity field in the QHD equations we need to present the many-particle 
wave function via two real functions as 

i;(R,t) =a(R,t)exp( iS ^ ). (7) 



h 

Average of square of the amplitude of the wave function a 2 gives the particle concentration, 
then the gradient of the phase of the wave function on coordinate of i-th particle gives 
velocity of i-th particle 

Vi(iM) = —ViS(R,t), (8) 
rrii 

which actually depends on coordinate of all particles in the system. 

Traditionally hydrodynamics equations contain velocity field v(r, t), which is velocity of 
the local center of mass, instead of the momentum density. These quantity simply relate 

j(r,t) =n(r,t)v(r,t). (9) 

Having the velocity field and the velocities of each particle we can introduce the difference 
between velocity of each particle and the velocity field u,;(r,i?, t) = Vi(R,t) — v(r, t), what 
gives the thermal velocities. 

Putting the presentation (JTj) of the wave function into the momentum flux tensor defini- 
tion ([6]) and expressing result via both the velocity field and the thermal velocities we get 
the momentum flux tensor in the following form 

where p a/3 is the thermal pressure, and T a/3 is the quantum Bohm potential, which is the 
quantum pressure caused by the quantum nature of particles as the de-Broule wave. 
Explicit form of the thermal pressure is 

N 

p a P(r,t)= dRj2 S ( r - r i) a2 ( R ^) m iUfu^, (11) 
J i=i 
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as we may expect the thermal pressure depends on the square of the thermal velocities. 
General explicit form of the quantum Bohm potential arises as 

ft 2 f , „ x o / „ x d 2 In a 



m J i=i 



dx ai dxpi 



(12) 



This general formula does not allow us to connect it with the particles concentration and 
the velocity field to get closed set of the QHD equations. In order to obtain the closed set 
of the QHD equation we consider this term in approximation of independent particles. In 
this case we find that the quantum Bohm potential in the following form 

T afj (T,t) = -^-(<PdPn(r,t) - -i-^(^n(r,t))(^n(r,t))Y (13) 

where we find that the quantum Bohm potential expressed via the spatial derivatives of the 
particle concentration. Divergence of this tensor has well-known in literature form 



Now we have represented the hydrodynamical quantities via velocity field and got explicit 
form of the quantum Bohm potential we can rewrite the set of QHD equation 

$n(r,t) + V(n(r,t)v(r,t)) = 0, (15) 

and 



inn 



(r, t){d t + vVK(r, t) + d^(r, t) + T a/3 (r, t)) = F a (r, t). (16) 



This equations have traditional form. Now we should consider the right-hand side of the 
Euler equation (I25j) . where we have the force field density. We have deal with the system 
of neutral particles interacting by means of the short-range interaction potential. We can 
use this main properties of the interaction to construct a general hydrodynamical theory for 
such system. 

III. SHORT RANGE INTERACTION IN QHD 

Let's start studying of the explicit form of the force density. We have not presented 
explicit form of the force density obtained in equation (jHJ). The force density consists of 
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two parts, which are action of an external field and an interparticle interaction, thus we can 
write 

F = -nW + F mf (17) 
and the second term, which is under our main interest, has following form 

F mt = - J di?^ 5(r-r i )(V i C/(r ii ))^ + (i?,t)^(i?,0- (18) 
There is no difference between i-th particle and j-th particle, thus we can rewrite it as 

F int = - J dR^SiT-TjiVjUiTi^+fatMl^t). (19) 

Using symmetry of the wave function to permutation of arguments and Tj, and the fact 
that Vjf/(r,j) = — ViU{vij) we can rewrite formula f fl9|) as 

F mt = f dRj28(r-r3)(ViU(rij))4> + (R,t)4>(R,t). (20) 

Keeping in mind that formulas ffT8j) and ( 120]) are different forms of the same quantity Fj nt . 
Thus, we can present F int as half of the sum of ({TBI and (|20|) . and consequently we find 



F 



int 



~\ / dRj^(5(v-r l )-5(v-v J ))^ l U(r, l3 ))ij + (R,t)^(R,t), (21) 

This formula is very useful for further transformation and let us proceed in (l2Tj) to variables 
of the center of gravity and variables of the relative distance of the particles defined as 

Rij = i(rj + r,-), Tij = Ti- rj ■ (22) 
Putting new variables in the force density F int we obtain 



F 



in t 



\ J dR E ( 6 ( r - - - s ( r - R ^ + \ t *)) (yvM) 



^ + (i?Ar_2, Rjj + -r^Ry — -Tij,t)^j(RN-2, Rij + 2 r *i'^i — T^ 1 "*- 7 '^' (23) 
where R^- + -r^ stands for i-th argument of the wave function, R^ — |rjj is the same for 
j-th argument, and Rn-2 is the set of rest coordinates. 

Since the interaction forces between the particles rapidly descend at distances of the order 
of the interaction radius, small \rfj\ give the main contribution to integral in (121]) . There- 
fore, in expression fT2Tl) . we can replace the multipliers at the interaction potential by their 



expansion in series by r^. It is worthwhile to admit that we should make expansion of three 
functions, they are the difference between 5 functions, i-th and j-th arguments in ip*(R,t) 
and ijj(R,t). We have reached the conclusion that the density of the interaction force for 
bosons with a short-range interaction potential can be represented in the form of divergence 
of the tensor field <9^cr a ^(r, t). Here, a a ^(r, t) is the quantum stress tensor conditioned by 
the occurrence of interparticle interaction. Divergence of this tensor is represented by the 
formula 



^4K( r -R, 



*j *j u 



x + \r% (V K + - <9£_)V + c.c. 

+ (Vk + <9£ + - 2d^_ + d^Ji, + c.c.) 

+ (^ + -^J^ + (^r + -5rJ^, (24) 

where Ry stands for z-th and j-th arguments, 9r + and <9r_ are derivatives on i-th and j-th 
arguments correspondingly. In this expansion we have included three first term of the series. 

This expansion is valid due to the fact that an interaction radius is small. The first 
(second, third) term of this expansion is proportional to the interparticle distance limited 
by the interaction radius in the first (second, third) degree. Therefore we can call this 
expansion as expansion in series on the interaction radius. Consequently, the first term in 
formula f l24l) arises in the first order by the interaction radius (FOIR). 

n 

We follow the basic paper [18(, where the QHD method was developed up to the third 
order by the interaction radius (TOIR). 

Using notion of the quantum stress tensor we can rewrite the Euler equation ffl6l) as 

mn(r,t)(d t + vV)v a (r,t) + dJ p a/3 (r, t) + a a/3 (r, t) + T a/3 (r, t) ] = -— n(r, t) V a V ext (r). (25) 

\ J m 



221 ] giving 



In the classic hydrodynamics of neutral particles a generalized Euler equation 
the Navier-Stokes equation at expansion of the stress tensor in the Taylor series on the 
spatial derivatives of the velocity field, restricted by linear terms only has form 

mnd t v a + (vV)/ + d p P aP = 0. (26) 
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In contrast to the general classic Euler equation our derivation gives explicit inner structure 
of P a/3 , that leads to more careful understanding of different physical mechanisms contribu- 
tion. 



IV. QUANTUM HYDRODYNAMICS OF BOSONS IN THE FIRST ORDER BY 
THE INTERACTION RADIUS 



We have obtained general form of the quantum stress tensor presented by formula 
In this section we expand further transformation of the first term in formula (1241) . Let's 
start transformations with rewriting of the first term in expansion of the quantum stress 
tensor. 



-,*) = -\l dR Y. *(r - R,,) | S i f { {rij U + (R,W(R,t). (27) 

J ..... Mi O Fjj 

In the integral (127|) dR can be rewritten as dR^-2dvidvj explicitly distinguish integration 
on coordinates of 2-th and j-th particles. Using variables of the center of gravity and the 
relative distance of i-th. and j-th particles we write dRN- 2 dridrj = dRN- 2 drijdIi,ij. It is 
also important to admit that both i-th and j-th arguments of the wave function in formula 
(12 7p are equal to Rij, after the expansion in a series on r^. Now we can see that integration 
on and R^ separates, and we can rewrite formula (1271) in the following form 

o*(T,t) = ~Tr(n 2 (r,r',t)) J d/-^^^, (28) 

where 

Tr/(r,r') = /(r,r), 

and we have also used notion of the two-particle concentration, which general definition is 
n 2 (r,r',t) = J dR ^ S(r - n)6(T> - r,)^*(i?, t)tp(R, t). (29) 

This definition can be rewritten in more useful form 

n 2 {r,r',t) = N(N - 1) J dR N - 2 {ni,n 2 , ... \r,r' ,R N _ 2 ,t){r, r',R N _ 2 , t\ni,n 2 , . . .), (30) 
where (ni, n 2 , . . . |r, r', Rn- 2 , t) is the N-particle wave function in representation of the oc- 

N 

cupation numbers, and dRx_ 2 = Y\_ dr k- 



a Q/3 (r, 



fc=3 
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For further transformation we need to extract evolution of particles related to arguments r 
and r'. For this purpose, we consider expansion of the wave function (r, r', Rn-2, t\ni, n 2 . . .) 



23l | . In the case of bosons, making expansion on one of arguments, we find 

(r,r',R N _ 2 ,t\n 1 ,n 2 . . .) = ^ J (r,t\f) (r', R N _ 2 , t\m, . . . (n } - 1),. . .), 

where we have that particle in an arbitrary quantum state / gives dependence on the coordi- 
nate r, and all particles alternately make contribution in (r, r', Rn-2, t\rii,n2 ■ ■ •) via (r, t\f) 
due to summation on all states. 

Making the second expansion of the wave function, including symmetry of the wave 
function due to permutation of arguments, we obtain 

(r,r',R N - 2 ,t\n 1 ,n 2 . . .) 

= £ E y^A/^^^tl^^tl/Ox^^tK,...^-!),...^-!),...) 



+ E y < r ' < r '' * ^ (^-a. ■ . - (n/ - 2) , . . .>. (31) 

where (r, t\f) = <ff(r,t) are the single-particle wave functions. Formula (T5TT) consists of 
two terms. The first term describes contribution of two particles from different quantum 
states. The second term gives contribution of two particles being in one quantum states, 
that corresponds to the fact that several Bose particles can exist in a one quantum state. 

Calculation of the two-particle concentration requires integration of the product of two 
wave functions in formula (1301) 



(m, . . . (ny - 1), . . . (n f - 1), . . . \n h . . . {n q > - 1), . . . (n g - 1), . . .) 

= *(/ - q)8(f ~ 90 + 5 (f ~ l'W - g), (32) 

and 

(m, . . . (n f - 2), . . . |m, . . . K - 2), . . .) = 5(f - q). (33) 

Formula (132|) explicitly reveals symmetry of bosons wave function. The second term gives 
exchange term, and, consequently, contribution of the exchange interaction. Using these 
formulas, after some calculations, we find following result for the two particle concentration 

n 2 (r, r', t) = n(r, t)n{r', t) + |p(r, r', t) | 2 + ^ n g {n g - l)\(p g (r, t) | 2 |^(r', t) | 2 , (34) 

9 
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where n g is a number of particles in the quantum state <p g , with a set of quantum numbers 
9, 

n(r, t) = Y^ %¥^( r ; t)(p g (r, t) (35) 

9 

is the particle concentration in terms of the arbitrary single-particle wave functions (p g (r,t), 

p(r,r',t) = ^n^(r,t)^(r',t) (36) 

9 

is the macroscopic density matrix. 

Putting (pUJ) in formula (1251) for the quantum stress tensor we obtain 

c^(r, t) = ~T8^ ^2n 2 (r, f) + p B (v, t)\ , (37) 

where 

p B (r,t) = ^rigirig - l)(|<^(r,t)| 2 ) 2 , 



and 



T = ^/*(r,.^) (38) 



is the interaction constant appearing from the integral on the relative distance of the parti- 
cles. Formula (135|) arises for the spherically symmetric interaction potential. In this chapter 
we have been developing the theory for neutral particles in the Bose-Einstein condensate 
state. One of the powerful and famous methods of this studying is the Gross-Pitaevskii 
equation, which we will derive below from presenting method. Integrating by parts in for- 
mula (1351) and assuming that the potential satisfies the condition that the quantity r 3 U (r) 
tends to zero as r tends to zero and infinity, we obtain 

T = -J drU(r), (39) 

that corresponds to the Gross-Pitaevskii result. 

To investigate solitons in the BEC, we use the set of the QHD equations up to the TOIR 
approximation [l^. The calculation of the first member in a quantum stress tensor that 
corresponds to the GP equation is made in [l8| under the condition that the particles do not 
interact. A more complete investigation into the conditions of the GP equation derivation 
from the MPSE shows that the GP equation appears in the first order by the interaction 
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radius (FOIR), if the particles are in an arbitrary state that can be simulated by a single- 
particle wave function. Such a state can particularly appears as a result of strong interaction 
between the particles that takes place in the quantum fluids. 



p B (r,t) = ^n g (n g - l)|^(r,t)| 4 

9 

= E^I^( r ' t )| 4 = iV2 |^o(r,t)| 4 

9 

= (iV|^ () (r,t)| 2 ) 2 = n|(r,t), (40) 
Finally, the quantum stress tensor in the BEC state has form 

af EC (^t) = - 1 -T5^nl^t). (41) 

We can see that the quantum stress tensor depends on the constant of interaction T and 
square of the particles concentration. 

We have been calculating the quantum stress tensor, which describes interparticle inter- 
action in the Euler equation. We have found cr^ EC , thus we can put in the Euler equation. 
Superfluid motion is the eddy-free motion. This means that v = Vip. For the isotropic ki- 
netic pressure p a/3 = p6 a/3 under barotropicity condition we introduce the chemical potential 
fi(r, t) as 

VMM) = 25fc0 (42) 
mn{r, t) 

This designation corresponds to the Gross-Pitaevskii equation. However, in the thermody- 
namics and hydrodynamics Vp(r, t)/mn(r, t) is equal to the gradient of the enthalpy, at the 
barotropicity condition. 

In the result we come to the Euler equation in following form 

md t v a {r, t) + -m<9V(r, t) + m<9>(r, t) 



- — d a ^h SA - Td a n(r,t) = -d a V ext {v,t) (43) 

We have obtained final form of the Euler equation. At description of the BEC dynamics it 
is enough to use the couple of the continuity (fT5|) and Euler (143!) equations, which have been 
obtained. As it was mentioned the Gross-Pitaevskii equation is usually used for the BEC 
description. The Gross-Pitaevskii equation is the non-linear Schrodinger equation and it 
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might be easily transformed in the couple of equations having hydrodynamical form. These 
equations coincide with equations obtained above and now we have the inverse problem. We 
need to derive the non-linear Schrodinger equation from the set of the QHD equations to 
present consisted description of the BEC. 

As the first step of the non-linear Schrodinger equation derivation we show existence of 
the Cauchy-Lagrangian integral for the Euler equation (1431) . Using notion of the potential of 
the velocity field we can rewrite equation ( |4"3]) as the gradient of the scalar function, which 
is the Cauchy-Lagrangian integral 

fl^(r, t) + ^(V0) 2 (r, t) + Mr, t) - -Tn(r, t) - ^1 A V^M) + (r t) = const. (44) 

Equation (144j) is the equation of the potential of velocity field evolution. 

Having equations for the particle concentration and velocity field potential evolution we 
can derive equation for function $(r, t), defined as 

$(r, t) = ^n(r,t)exp (^m0(r, t^j , (45) 

which called the macroscopic wave function, the order parameter or the wave function in 
the medium. This construction has macroscopical meaning due to it's definition via the 
macroscopic parameters. 

Differentiating the macroscopic wave function ( 1431) with respect to time and using the 
continuity equation (fT5"j) and the Cauchy-Lagrangian integral (|44p we find following non- 
linear Schrodinger equation 

thdtQfat) = (- T ^ + fi(r,t) + V ext (r,t)-T\ $(r,t) | 2 )$(r,t), (46) 

which is well known as the Gross-Pit aevskii equation 

a, a. 

Operator ft 2 V 2 /2m arises 

in the Gross-Pitaevskii equation due to the quantum Bohm potential ffl2l in the Euler 
equation (1251) . or more precisely, it appears due to the approximate form of the quantum 
Bohm potential. 

The wave function $(r, t) is normalized by the condition 

J dr$(r,t)*$(r,t) = N, 

where N is the number of particles in the system. 
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V. CONTRIBUTION OF THE TEMPERATURE IN DYNAMICS OF BOSE PAR- 
TICLES IN THE FIRST ORDER BY THE INTERACTION RADIUS 



We describe QHD description of bosons at non-zero temperatures. Other method of the 
BEC description, for example the self-consistent Hartree-Fock-Bogoliubov approximation, 
can be found in Ref.s [scj-jssj]. 

At low, but non-zero, temperatures part of Bose particles are in the BEC state (on the 
ground energy level) and other part distributes on excited states (non-condensed particles). 
Traditionally one species of bosons in described conditions is considered as a mixture of two 
liquids. One liquid is the particles in the BEC state, and another one is non-condensed 
particles. Since we consider non ideal Bose gas, we have interaction between particles, and 
in chosen model we have as interaction between particles of each liquid and inter-liquid 
interaction. 

Interparticle interaction can lead to exchange of particles between the two liquids (since 
these liquids contain atoms of the same species), even if we keep system at fixed temperature. 
Decreasing (increasing) of the system temperature leads to decrease (increase) of a particles 
number in excited states, and, correspondingly, it leads to increase (decrease) of a particles 
number in the BEC state. In mentioned cases the number of particles in each liquids changes, 
but total number of particles does not change. Therefore we have the continuity equation 
for total particles concentration n, which has usual form 

«9 t n(r,t) + Vj = 0, (47) 

where j is the current of all particles, or the total current, but for described above processes 
we have creation and destruction particles, at transition between liquids, which cause addi- 
tional terms in the right-hand side of the continuity equation corresponding to the partial 
concentration introduced for each liquid. 

We will consider a system of Bose particles at fixed non-zero temperature and neglect 
by the particles transitions between the liquids. In this case we the continuity equation for 
each liquid in usual form 

d t n B (r, t) + V(n B (r, t)v B (r, t)) = (48) 
for particles in the BEC state, and 

d t n n (r, t) + V(n„(r, t)v n (r, t)) = 0. (49) 
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for the non-condensed bosons, where the total particles concentration n(r, t) and current 
j(r, t) are divided on two parts 

n(r,t) = n B (r,t) +n n (r,t), 

and 

= j B (r,t) + j n (r,t), 

ris, n n are particle concentrations for the BEC and non-condensed particles, and j B and j n 
are corresponding currents, which have usual relation with the corresponding velocity fields 
j B = n BV B and j n = n n w n . 

Next step is dividing of the momentum balance equation (the Euler equation) on two 
parts corresponding to condensed and non-condensed particle evolution. 

We start with kinematic part of the momentum evolution which is given by the momentum 
current IT a/3 . It has bilinear structure on wave function as the concentration n and particle 
current j. Thus we can represent the momentum current as the sum 
Explicit form of II^, and TYjf are given by formulas ([TO]) and f|T2|) , where we should put n b , 
v B or n n , v n instead of n and v, and neglect thermal pressure p a/3 . 

For understanding of dynamical part of momentum evolution evolution of BEC j#(r, t) 
and non-condensed particles j n (r,t) we need to consider formulas ( 13T1) and (1341) in details. 

The first multiplier in formula (13TI) . which has argument (r, t) is related to the particle 
whose motion we consider. Other one particle wave functions are related to the particles that 
influence on dynamic of considered current. This is give us ability to obtain the separate 
equation of dynamic atoms in the BEC state and the non-condensed state. If one-particle 
wave function with argument (r, t) describe the BEC state (has subindex "B"), we put this 
term in momentum balance equation for the BEC. In the case one-particle wave function 
with argument (r, t) describe the non-condensed state we put this term in the momentum 
balance equation for the non-condensed particles. 

If we consider dynamic of particle in the BEC state, it means a quantum state / in 
formula (T5T|) describes the BEC state we find that the first two terms in formula (j34p 
describes quantities related to the particles in different state, for our case it means that 
they describe interaction of condensed and non-condensed particles. Thus the square of 
concentration and the square of module of the density matrix lead to evolution of the BEC 
due to interaction with non-condensed particles, they trace arises as 2n n n B . The last term 
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in formula (}3Tj) describes two-particle in the same quantum state /. It gives us that the 
last term describes interaction between particles in the BEC state, and as we obtained in 
formula (HUj) it equals to the square of the concentration of condensed particles. 

Let's consider momentum evolution of non-condensed particles. Now a quantum state 
/ describes one of the non-ground states. We repeat that the first term (two first terms) 
in formula (l3Tj) (in formula (|34j) ) describes quantities related to the particles in different 
quantum states. In the case of the BEC evolution we had fixed state of evaluating particles. 
Thus we had only one combination of of different states they are BEC state and some of non- 
condensed states. Now we have more possibilities. Considering particle in a non-condensed 
state we find that it can interact with a particle either in the BEC state or in other non- 
condensed state. Thus we obtain the sum of 2n n riB and 2n l n . The last term gives additional 
contribution in non-condensed particles interaction. However, it was shown in Ref. 26| that 
it gives small contribution and we can neglect this term. 

As a result we obtain 

h 2 

mn B {d t v% + viv p v%) d a An B 

Am 

+ ^ { ^ nB n ^ nB ) = -n B V a V ext (r, t) + \rd«{2n B n n + n|), (50) 



and 



h 2 

mn n (d t v a n + ^V^<) - — d a /\n r , 

Am 



+ ^{ p^ d " Un ) = -n n V a V ext (v,t) + \r&*{2n B n n + 2rfl, (51) 

These equations have similar form, especially due to kinematic nature of the left-hand sides 
of the equations. The right-hand sides describe interaction with external field V ex t and 
interparticle interaction proportional to the interaction constant T. We have only one inter- 
action constant since we consider one species in different quantum states. The right-hand 
side of equation fl50l contains two terms proportional to T. One of this terms contains 
n 2 B and describes interaction between particles in the BEC state. The term, containing 
n B n n , describes interaction between non-condensed particles and particles in the BEC state. 
Analogous terms in equation (15TT) have similar meaning, and term n B n n present interaction 
between two fluids, and it's influence of non-condensed particles evolution. The last term in 
equation fl5TT) . proportional to 2n 2 n , describes interaction in non-condensed particles. 
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The last term in equation (150]) emerges from the last term in formula (I34p . and the last 
term in equation (j5~Tj) appears from the two first terms in formula fl3l|) . therefore it is quite 
normal that they have different coefficients. 

VI. QUANTUM HYDRODYNAMICS OF BOSONS AT ACCOUNT OF INTER- 
ACTION UP TO THE TOIR 

We have expanded derivation of both the Euler equation and corresponding non-linear 
Schrodinger equation, which is the Gross-Pitaevskii equation. In this section we will briefly 
discuss a contribution of the second and third terms in the quantum stress tensor. 

In the second and third terms, as in the first one, the variables of the center of gravity 
and variables of the relative distance of the particles can be separated. The second term in 
the stress tensor, for spherically symmetric interparticle potential, equals to zero due to the 
integral on the variable of the relative distance. Thus, we need to consider the third term 
only. 

In the absence of particles in excited states, the quantum stress tensor can be presented 

as 

<c(M) = -l'T8* > 52n g (n 9 - l)(|^(r, *)| 2 ) 2 

g 

~T 2 (S^A + 2d a d?) - l)(|^(r,t)| 2 ) 2 

g 

-T 2 (rt 75 + 5^5^ + 5 aS 5^)x 
x ^ n g( n g ~ !) ^P*g ( P* g ( < Pg d 'Y d 6<P9 - 9j<fgd s <fg) + h.c^j , (52) 
where we assumed that particles occupy a ground quantum state described by wave function 

<Pg- 

The quantum stress tensor in the FOIR approximation was found for system of particles 
which are in some quantum state described with wave function (po(r,t) (see formula ( l40l) ). 
It was managed to calculate the quantum stress tensor, using intermediate formula ( 1521) . in 
the TOIR approximation for approximately independent particles, as for free particles and 
for particles in the parabolic trap [3]. In this case the second term in formula fl52l) simplifies 
and the last becomes equal to zero. 
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Method of calculation is the same as described above for the first term, but rather more 
complicated, so we do not present details. Final form of the third term in the quantum 
stress tensor is 



where 



a Q ^(r, t) = -^T 2 (5 a/3 A + 2d a d?)n 2 (r, t), (53) 



T 2 ee ^ I (rf^-dr (54) 
30 y dr K J 



is the constant of interaction arising at account of the short-range interaction up to the 
TOIR, this definition differs from the definition presented in Ref. 
1/8 in the integral ([51 



18j , here we put factor 



Collecting parts of the quantum stress tensor arising in the FOIR and TOIR we find 

a^(r, t) = --TcTV(r, t) - -T 2 (5 a ^A + 2d a d^)n 2 (v, t), (55) 
2 6 

which has form of an operator acting on the square of the particle concentration. Formula 
( )55l) reveals symmetry of the quantum stress tensor: a a P = a^ a , which exists in general 
formula ( )2~4"1) and remains in approximate formula ( 1551) . 

Euler equation contains divergence of the quantum stress tensor which emerges as 

d^(r, t) = -Tn(r, t)d a n(r, t) - ^T 2 d a An 2 (r, t). (56) 

We have found that the quantum stress tensor depends on higher than first spatial deriva- 
tives, see the second term in formula (|56|) . Such dependence leads to nonlocal non-linear 
Schrodinger equation. Other nonlocal approximations for the BEC was considered in Ref.s 



2i, |25|. 

At studying of quantum gases we have dial with boson-fermion and fermion-fermion 
mixtures. The QHD of ultracold fermions and boson-fermion mixtures was developed up to 
the TOIR approximation in Ref. 181 ] . 



VII. ENERGY EVOLUTION UP TO TOIR APPROXIMATION 

Energy evolution is associated with both the collective motion and thermal motion of 
particles. We have no thermal motion in the BEC. Thus, it was enough to use two hydro- 
dynamical equations for the BEC description. Non-condensed particles are involved in the 
thermal motion. At low enough temperatures we can neglect by the thermal pressure in the 
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Euler equation for non-condensed particles. Including thermal pressure contribution and 
influence of it's evolution we have to have equation for the pressure evolution. In most cases 
we can limit our treatment with consideration of the scalar pressure which relates to the ki- 
netic energy of the thermal motion. We actually can derive the pressure evolution equation, 
but it will be more complicated due to tensor nature of the pressure, so we consider energy 
evolution only. We give definition of whole energy, including energy of collective motion 
described fully by the momentum density evolution, and we will extract energy of thermal 
motion below. 

Energy density for quantum system is defined as 

e(r,f)= / dR^6(T-T i )(±L*pfy + c.c?) + \f^U(\T i -T S \)w\ (57) 

i=l \ 1 ' J 

where the first term in the big brackets is the kinetic energy density of i-th particle, and the 
second term is the density of potential energy. 

Differentiating the energy density (|57|) with respect to time we find general form of the 
energy evolution equation 

Qte(r,t)+VQ(r,t) = -j(r, t)VV ext (r, t) + A(r, t), (58) 

which has traditional structure and does not distinguish from classical case. 

In equation (15 8j) we have the energy flux Q(r, t) and the density of the work A(r, t). Due 
to our calculation these quantities arise in explicit form. Obtained equation is appropriate 
even for particles with long-range interparticle interaction, so explicit form of Q(r, t) and 
A(r, t) are correct in general case. 

Let us present explicit form of the energy flux Q(r, t), which we separate on two parts 
having different meaning 

Q a (v,t) = Q? kin) (T,t) + Q% nt) (v,t), (59) 

the flux of the kinetic energy Q" fcin )(r, t), and the flux of the potential energy Qf int )(r, t). 

The flux of the kinetic energy has following representation via many-particle wave- 
function 

Q?kin)(*,t) = I dR^6(r-r t )^((ptijrp^ + rPt& + cxX (60) 
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the flux of the potential energy containing the interaction potential Uy arises as 

Q? mt) (r,t)= 1 - J di?^5(r-r i )^^UW + (P^)*A (61) 



and the work density is 



i r N ( 1 

A(r,t) = -- I dR^SiT-TiXVfuJ — 



1 



Using the fact that we consider short-range interaction, thus we can use the method 
described in section (III) to represent the work density and the potential energy flux in 
corresponding form. This form appears to have very large form, so we present these quantity 
separated on parts. We start with representation of the potential energy flux via one-particle 
functions ip g describing one-particle states, which can be occupied by particles. Dividing 
Q^ nf) on two parts 

Q%nt)( Y ^) = Q(int)d( r ^) + Q(int)s ( r > *) " 

Indexes d and s mean that particles under consideration are in different states and same 
state, correspondingly. 

Explicit form of Q? int \ d arises as 

Q(int)d — ( E n 9 n 9 ,< P*g < P*A da 'P9) < Pg' ~ C - C - J 

\g,g',g^g' / 

+ i r2 ( 2 E n g n gl ^ gl {^d^ g )<p gl -2 n 9 n g ,^;,(d a d^ 9 )d^ 9 , 
V 9,g',9¥=g' g>g',g^g' 

- E w((9>p^(^v,)^+^(9v*o(^j^ 

9,9',g¥=g' V 



+2 Yl n g n g/ cp* g (d a cp* gl )(d^ g )d^ gl - c.c. 

9,9',9^9' 

+ jLr 2 a/W y 

2im \ 

\9,9 ' ,9¥=9' 



g,g',g^g' 

n 9 n g><P*a<P* g >( da( Pg) l Pg' ~ c - c - I ( 63 ) 



\g,g\g+g' 

here and in the following formula for Q^ int ^ s we use following designations 

Fi = / dxU(r), (64) 
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and 



where 



rf = J rVf/(r)dr = 5 a/3 T 2 , (65) 
Air f 

T 2 = — U{r)r 4 dr. (66) 

Coefficient T\ describing interaction in the FOIR approximation is the Gross-Pitaevskii 
interaction constant g, which is simply related with T = —g = —V. T 2 arises in terms 
coming out in the TOIR approximation, so it has to be connected with the T2. Integrating 
definition of T2 (|54l) by parts we find that T2 = — 1^. 

The potential energy flux related to interaction of particles in the same quantum state 
has form of 

Q(int)s(r,t) = [^n g {n g - l)<p*<p*<p g 8Ptp g - c.c. j 

V 9 9 

"E^K - l)(d>>>^V + - l)<P* g (d a <P* g )(d^ g )d^ a - c.c. 

9 9 1 

+ (E %K " 1 W r . ^> 9 ^ 9 ~ cc^j ■ (67) 

We are interested in studying of particle dynamics at low temperatures, consequently we 
suppose that the ground energy state and an interval of low energy states are macroscopically 
occupied. We also can approximately neglect tail of states with larger energies which are 
faintly occupied. Therefore we will consider cases then n g {n g — 1) « n 2 g in Qj nt ( s ) and A s . 

We have explicitly presented intermediate formulas for the potential energy flux up to 
the TOIR approximation. Now we consider analogous formulas for the work. 

We divide the work on two parts following the recipe as for potential energy flow 

A(r,t) =A d (r,t)+A s (r,t), 

A s (Ad) is the work related to two in different (same) quantum states. 
Following formula for the work A^ emerges from our calculations 



M T i t) = ^ T9a [ E n 9 n 9'^g'^9 lda ^9 ~ C - C - 
y9,9',9¥=9' 
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+-^-T 2 (5 Q/ V' 5 + 5^8^ + 6 aS 5^)d^dW 5 ( E n g n g ,<p*<p* g ,<p g ,d a <p g - c.c. 



^g,g',g^g' 



+-^-T 2 (rt^ + 8^8 ps + 5 aS 5^)dP ( E n g n g ,<p* g <p* g ,<p g ,d a dr'd 5 <p g 



^9,g',g¥=g' 



E Vs'^^)^!'- 2 E n 9 n 9' ( P*g ( P*A d5< Ps') aPiare 'Pa 

g,g',g¥=g' g,g',g¥=g' 



g,g',g¥=g' 

+ E w((^P^(^ 9 05V 9 + (9>*0^(5 7 ^)^V 9 ) -c.c). (68) 

9,9' ,97V V J ' 

In terms proportional Y 2 we meet the tensor structure 5 a/3 5 7<5 + 5 a7 5 /35 + (5 Q!<5 5 /37 , which has 
obtained in the quantum stress tensor cr a/3 (152j) . in the terms which appeared in the TOIR 
approximation. It relates to the fact that the second interaction constant T 2 emerges as 
a fourth rank tensor, and the structure under discussion reflects symmetries of this tensor. 
We also meet this structure in following formula for A s , which is 



We 

V 9 



M r >*) = 4^ T( M E^K - 1 )<P*g t P*g<Pg 



+^7^-T 2 (rt 7 < 5 + 5^5? s + 6 a5 6^)d^d^d s 1 E^K ~ l )v* a ^gd a Vg ~ c.c. 



2 3km 

n 



V 9 

+ E^K - i)<p m B <pl{8rt&<p g )ff><p g - 2 E%(% - i)<p;<p;(&< Pg )ff*8r'<p a 

9 9 

- 2E%K - lWPipDtpgffr&ipg + 2E%K - l)<p;(ff'<p;)(d'<p a )tf<p 8 - c.c. j . (69) 

9 9 / 

These expressions ( l63l) . ( |67j) . ( |68l) and ( |69l) can be used for developing of different approx- 
imations for the quantum gases evolution, and we use them for one simplest approximation 
of approximately free particles. Thus, we use plane waves 

<f P ( r , t) = -j= exp ^-^(M ~ P r )) ' ( 70 ) 
for description of one-particle states. 
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To find closed description for the energy evolution we need to present the particles current 
j a , momentum current U a/3 

r (r> *) = ^E n « t)p a ^ *) + c - c ) > ( 71 ) 

IT*(r, t) = ±- n 9 L* g (r, t)f/%(r, t) + (p a ^* g (r, t))fcp g (v, t) + cc.) , (72) 

£ = 4^ E n 9 K( r ' ^V 9 (r, t) + cc. J , (73) 

g ' 

Q a = ^E n /(^)*(r,t)pV,(r,t) + <p;(T,t)f?<p g (T,t) + c.c.Y (74) 



and 



These quantities get following form for the plane wave function (ITUl) n = XpTip/V, 
j a = X p p a n p /V, e kin = S pP 2 n p /y, = Z p p<*pPn p /V, and Q Q = E p p a p 2 n p /V. These 
formulas allow us obtain following formulas for the potential energy flux and the work 

Qf int)d = -Tnf - -T 2 A(nj a ) - -^{me {kin )f - nQ a ), (75) 

Q(int)s(v,t) = -±-TA a (r,t) - ^T 2 AA a (r,t), (76) 
A d = Td a {nf) + T 2 9 a A(nj a ) + ^Ild a (fe km + /IT* - 3nQ a km ), (77) 

and 

A s (r,t) = -^T9 a A Q (r,t) + — — <9 a AA Q (r, £), (78) 
2m 2m 

aim) = ^a;. (79) 
p 

Fixing local macroscopic quantum state via knowledge of macroscopic quantum param- 
eters n, j, £ we limit a set of microscopic quantum states, which needs for description of 
particle system. Overwise, if we know microscopic quantum state of a system (we know in 
which quantum state each particle is) we uniquely obtain macroscopic parameters describing 
system. 



where 
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VIII. THREE PARTICLE INTERACTION BETWEEN PARTICLES IN THE BEC 



QHD of the BEC with three-particle interaction was developed in Ref. [26|. Before 
presenting of the basic equation let us briefly discuss meaning of the three-particle interaction 
(TPI). We can present interaction between three particles V^rj, r^, r k ) in the form 

V 3 (v i ,r j ,v k ) = U(r i:j ) 

+ U(r lk ) + U(v jk ) + U{ r ik: r jk ), (80) 

where = rj — r^, U^ — U(\ |) is the binary interaction potential, Uj k — U(\ r^- |, | | 
, | r k j |) is the TPI potential that does not contain combination of binary potentials U^. It 
means that if we have deal with the simultaneously collision of three particles we can expect 
that the two-particle potential is enough for description of this process. If particle has no 
inner structure, as electron, for example, they have potential of two-particle interaction only 
If particles have an inner structure, as atoms and molecules, we can expect that the TPI 
Uij k can give contribution in the collision. 

The Hamiltonian of the system under consideration has the form: 

i i 

where f>1 = —ih'Vi is the momentum operator of the i-th particle, mj-is the mass of the i-th 
particle. Let us consider the TPI nonequivalent to combination of the binary interaction 
described by Uij. Interaction of three or more particles at the same time by means of binary 
potential is described by U^ and have no connection with the TPI. 

md t r(r,t) + d^(r,t) = - J dr\V a U(r, r'))n 2 (r, r', t) 

- j dr' J dr"(d a U(r,r',r"))n 3 (r,r',r",t)-n(r,t)V a V ext (r,t). (82) 
we have introduced the three-particle concentration 

:r,r',r",t) = J dR £ 5(v - r,) 
x 8(r' - rj )5(r" - r k )i; + (R, t)^{R, t). (83) 



n 3 [ 
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It is necessary to consider variables of the center of gravity and variables of the relative 
distance of the particles for three particles 

R-ijfc 3(^1 ~^~^k)j ^ij ^jj ^ik ^k ' 

T jk = r j — r k = r ik ~ r ij- (84) 

Using these variables we can represent the force field in terms of the quantum stress tensor 

d tJ a (r,t) + -^(IT^M) + ^(r,t)) = --n(r,t)V a V ext (r)- (85) 

Limiting our consideration by the FOIR approximation for two- and three-particle inter- 
action we find following representation for the quantum stress tensor 



\j dR £ ( 



4% + r ik d k)U{r lv r ik , I r - - r tt |) ] 5(r - R ijfc )V + ^, (86) 



9 

where 

ijj(R, t) = Rijki •••) -Rijfc, £) 

for ^ function in the term describing the binary interaction, similarly in the term describing 
the three-particle interaction the psi-function has following structure 

1p{R, t) = lf)(..., Rijk, Rijki •••) Rijki •••) t). 

Separating variables of the center of gravity and variables of the relative distance of the 
particles in the quantum stress tensor we obtain 



a aP {v,t) = - l -Tr{n 2 {vy,t)) j dr- 



~Tr{n 3 {r, r',r",t))x 



r dr 



x J dr 12 dr 13 (r^ 2 -^U(r 12 ,r 13 ,\ r 12 - r 13 |) + rf 3 ^[/(r 12 , r 23 , | r 12 - r 13 |)j 



(87) 
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where 

Tr/(r,r') = /(r, r), 

and 

Tr/(r,r',r") = /(r,r,r). 
Calculation of the three-particle concentration requires further expansion of the many- 



particle wave function in comparison with formula ( l3Tj) 26 1. 



For particles to be found in the Bose condensation state contribution of Trn^ reduces to 
m B (r,t) = ^n g {n a ~ !)K - 2)K(M)| 6 = n%(r,t), (88) 

9 

which was calculated analogously to (140|) . 
In the issue we get the Euler equation 

md t v a (r, t) + -md a v 2 (r, t) + m<9 a /i(r, t) 

-^&*An + fd^ ( ■ 9 " n ) - Td"n(r, t) 
Am Am \ n J 

- 2 X a V(r, t)d?n(r, t) = -n(r, t) V a K*t(r, t). (89) 
where the interaction constant for the TPI arises in the tensor form 

1 f j_ j_ / / r l r ln , r 2 r 2 a , oK - r 2)( r l - r 2) 



6 y \ \ ri r 2 ri - r 2 



x U (ri, r 2 , y r? + rf + 2r x r 2 cos fi) I , (90) 

where f2 is the angle between r 1 and r 2 , ri and r 2 are modules of the vectors and r 2 , 
and <9i, 9 2 , c?3 are derivatives of function U on its arguments r\, r 2 , y^rf + r\ + 2rxr 2 cos f2 
correspondingly. We can see that = X^ a - 

Reducing our description to the scalar TPI we can rewrite the Euler equation as 

md t v a (r, t) H — md a v 2 (r, t) + m<9 a /i(r, t) 

4m 4m \ n y 

- X 9V(r,t) = -9 Q ^(r,t), (91) 
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where scalar TPI constant has following form 



X = -3 I 



x U (ri, r 2 , y rf + r| + 2rir 2 cos fi) J . (92) 

Including of the three-particle interaction in the FOIR approximation gives the one ad- 
ditional term in the momentum balance equation and it also gives contribution in the non- 
linear Schrodinger equation, the method of derivation described above for Bose system with 
the two-particle interaction, which reveals as a nonlinearities of fifth degree 

ihdMr, t) = {-~^ + Mr, t) + V ext (r, t) - T | $(r, t) | 2 - X | $(r, t) | 4 \ $(r, t). (93) 

That is well known Gross-Pitaevskii equation with the nonlinearity of the fifth degree 
5], which arises due to the TPI, in approximation of the scalar TPI. 



IX. DISPERSION OF LINEAR COLLECTIVE EXCITATIONS IN BEC 

In this chapter we have considered the two-particle interaction up to the TOIR and the 
three-particle interaction in the FOIR. In this section we consider contribution of these 
interactions in the dispersion of the collective excitations of the BEC. Let's write the set of 
the QHD equations including described interactions, which is 

d t n + V(nv) = 0, (94) 

and 

mnd t v a + -mnV a v 2 - ^-d a An + ^-d? ( ^ ' ^ 
2 4m 4m \ n J 

- l -Td a n 2 - X d a n 3 - ^T 2 d a An 2 = -nV a V ext . (95) 

The last three terms in the left-hand side of equation (195)) describe interparticle interaction, 
which are two-particle interaction in the FOIR approximation, three-particle interaction in 
the FOIR approximation, and two-particle interaction in the TOIR approximation, corre- 
spondingly. 

We consider the small perturbation of equilibrium state like 

n = n + Sn, v a = + v a , (96) 
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Substituting these relations into system of equations ( 194)) and ( 195]) and neglecting nonlin- 
ear terms, we obtain a system of linear homogeneous equations in partial derivatives with 
constant coefficients. Passing to the following representation for small perturbations Sf 

5f = f(cu, k)exp(-itut + zkr) (97) 

yields the homogeneous system of algebraic equations. The magnitude of concentration 
of the BEC is assumed to have a nonzero value. Expressing all the quantities entering 
the system of equations in terms of the concentration of BEC, we come to the dispersion 
equation for elementary excitations 

\ Am z ml \ m ml 

The first and third terms corresponds to the Bogoliubov spectrum, which also appears from 
the Gross-Pitaevskii equation. The first term arises from the linear part of the quantum 
Bohm potential. The third term exists due to the two-particle short range interaction in the 
first order by the interaction radius. Account of the two-particle interaction up to the TOIR 
gives the second term, which leads to dependence of coefficient at k A on the equilibrium 
concentration uq. The last term is caused by the TPI. 

Nonlinear shift of the frequency of collective excitations was calculated up to the TOIR 
approximation j^j, change of the form of the bright soliton in the BEC due to account of 



the short-range interaction up to the TOIR was obtained in Ref. 



Moreover, it was 



shown that considering of the short-range interaction up to the TOIR gives new solutions 



29| | . namely new soliton solution was found. This soliton exists in the repulsing BEC, and 
it reveals clS 3X6 cl of compression, so it was called bright-like soliton, when the bright soliton 
exists in the attractive BEC. At usual conditions the repulsing BEC reveals the dark soliton, 
which is an area of rarefication. 



X. QUANTUM HYDRODYNAMICS OF THE DIPOLAR BEC 

One of the more excited topics in BEC studies is the electrically polarized BEC (EPBEC), 
which is the BEC of polar molecules. Quantum gases of particles having electric dipole 
moment add two interesting fundamental properties. These are long-range interaction and 
anisotropy of the interparticle interaction. 
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The Gross-Pitaevskii equation with the cubic nonlinearities has been used for studying of 
the unpolarized BEC with two-particle interaction. Taking into account the three-particle 
interaction we come to the non-linear Schrodinger equation containing the quintic nonlinear- 
ity along with the cubic nonlinearity. Thus, we have generalization of the Gross-Pitaevskii 
equation. The cubic and quantic nonlinearities play role of the potential energy of interpar- 
ticle interaction in the non-linear Schrodinger equation. Therefore, if we want to generalized 
the Gross-Pitaevskii equation for description of the EPBEC we need to add the potential en- 
ergy of dipole-dipole interaction. In this way was obtained the generalized Gross-Pitaevskii 
equation for the EPBEC fl-Q 

ihd t $(r,t) = l--^— + fj t (r,t) + V ex t(T,t)+g \ $(r,t) | 2 

+ (f f ^^-f/ 1 $(r '' t} I 2 ) $ ( r '^ (") 
In equation (199]) following designation are used: $(r, t) is the macroscopic wave function, \x 
is the chemical potential, V ext is the potential of external field, g is the constant of short- 
range interaction, d is the dipole electric moment of single atom, m is the mass of particles 
and h is the Planck constant divided by 2ir. The last term in equation ( 1991) describes the 
dipole-dipole interaction. 



We derive the QHD equations for the EPBEC 



171 ] , (36| , and (37J from the many-particle 



Schrodinger equation with the following Hamiltonian 



(100) 

The first term here is the operator for kinetic energy. The second term represents the 
interaction between the dipole moment df and the external electrical field. The subsequent 
terms represent short-range Uij and dipole-dipole interactions between particles, respectively. 
The Green's function for dipole-dipole interaction is taken as Gf? = VfVf (l/r^). 
Usually people use following Hamiltonian of dipole-dipole interaction 

h m = y^-dyl (ioi) 
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which is enough for mechanical description of several particles motion, but it does not enough 
for construction of the field theory. Thus, in this paper for interaction of electric dipoles we 
use the following, more general, Hamiltonian: 



H dd = -d a d?- ■ d?4. 
r 



There is well-known identity 



_ = 6 6rrjr_ + 47^ 

r r 6 3 

so we can see the difference between usually using Hamiltonian and one's used in this 
paper. The corrections of our selection followed from the fact that the equations obtained 
in the paper coincide to the Maxwell equations. At derivation of the QHD equation in 
self-consistent approximation lead to field equation 

VE(r, t) = -4vrVP(r, t), (103) 

but if we used Hamiltonian fllOip we would obtain 

VE(r,i) = yVP(M), (104) 

instead of f TI03|) . 

Electric field arising in the QHD equations in the self-consistent field approximation at 
using of the Hamiltonian (jlOip satisfy to the two Maxwell equation (11031) and 

V x E = 0. (105) 

Scheme of derivation of the QHD equations is the same as described in previous sections. 
Thus, we have following set of equations, which contains additional equations in compare 
with the unpolarized systems described below. The continuity equation 

d t n + d a {nv a ) = (106) 

has the same form as usual. 

The momentum balance equation for the polarized BEC has the form 

h 2 _ . h 2 „ a ( d a n ■ d?n 



mn 



(dt + vV)/ - — d a An + — d 13 



Am Am \ n 
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= Tnd a n + P f) d a E p , (107) 

and contains one additional term. New quantity has appeared in the momentum balance 
equation, which is the electric dipole moment density 



P(r,t) = J dR^8{r - r^drt* (R,t)iP(R,t). 



(108) 



It is obviously very important quantity for description of polarized BEC evolution, so that 
we should consider their evolution and influence of the polarization evolution on dynamics 
of hydrodynamic variables (particle concentration n and velocity field v). When we consider 
system of fully polarized dipoles at zero temperature we can rewrite polarization P = dn. 
At assumption that particles dynamic do not lead to change of the dipoles direction we 
can conclude that changing of polarization caused by changing of particles position, so 
polarization increase in area where concentration increase and vice-versa. The generalized 
Gross-Pitaevskii equation ( 1991) corresponds to described case. Generally speaking we should 
include dipoles direction evolution, and consequently more detailed polarization evolution. 
Even in ferromagnetics, where magnetic moments are kept parallel due to strong exchange 
interaction, there are waves of dipoles direction (spin waves). Therefore, we need to derive 
equation of the polarization evolution. 

Considering two kind of interaction, the short-range interaction and dipole-dipole inter- 
action, the last one is an example of the long-range interaction, we have dial with differ- 
ent approximations for each interaction. Approximations for short-range interaction were 
minutely considered in the beginning of the chapter. We have used the self-consistent field 



approximation 17J for dipole-dipole interaction. Briefly, this approximation corresponds to 
approximate representation of two-particle function as a product of corresponding one par- 
ticle function. More precisely, we can explain it on example of two-particle concentration. 
The self-consistent field approximation corresponds to the first term in formula f )3^|) . which 
is the general representation of the two-particle concentration in terms of one-particle states 

<Pg- 

As usual we differentiate the polarization ( 11081) with respect to time and using Schrodinger 
equation with the Hamiltonian f llOOp the equation of polarization evolution appears as 

d t P a (r, t) + dPR a P(r, t) = 0, (109) 

R a/3 (r,t) is the current of polarization. 
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Using a self-consistent field approximation of the dipole-dipole interaction we obtain an 
equation for the polarization current R a/3 (r,t) evolution 

d t R aP + <9 7 f iTV + iTV - P a v?vA - -^-d AP a 
\ J Am 1 

+ * (Qjp^ + mn\ = Lrgf /a + iL f!fi^ ! (110) 

8m J \ n n J m \ J m n 

where the second term, containing three terms in large brackets, in the left-hand side is the 

convective part of the polarization current evolution. The third and fourth terms describe 

an analog of the quantum Bohm potential. We have two terms in the right-hand side 

describe influence of interaction on the polarization current evolution. The first term in 

the right-hand side presents contribution of the short-range interaction, where we have T 

as the signature. The very last term of the formula (lllOp describes interaction of dipoles 

with electric field and includes both external electrical field and a self-consistent field that 

created by particle dipoles. This term contain numerical constant a. 

Considering monochromatic collective excitations in the EPBEC 

5f = f(u, k)exp(-itut + zkr) (111) 

choosing equilibrium condition as n = uq, v = 0, P = Po& z , and R a/3 = 0, including that 
equilibrium external electric field is E = Eoe z , we find that the dispersion dependence for 
the collective excitation in the EPBEC can be expressed in the form of 

co 2 = —(— + 4vra^ - -Tn k 2 
2m \ 2m 2 



± \\ ( ^Tn k 2 + 4ira^-^J - 8nTk*P 2 j , (112) 

quantities with subindex are constants. This formula shows that in the EPBEC exists 
two waves, due to two signs in front of the square root, instead of one wave existing in an 
unpolarized BEC (|98|) considered in previous section. 

Comparing obtained dispersion dependence for EPBEC (11121) with the one derived in the 
absence of the electric dipole moment (I98jl for T2 = and x — 0, since we have not considered 
the TPI and two-particle interaction in the TOIR approximation, we derive Po — * limit of 
formula (I112p . which appears as 

2 1 fh 2 k A „ , 2 %1iP 2 k 2 \ , , 

oo 2 = - Tn k + °— , 113 

m \ Am n J 
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for minus in front of the square root in formula f II 1 2 j) and 

2 1 fh 2 k A 1 2 8vr(a - l)P 2 fc 2 



^ = ~ K— " ^Tnofc 2 + ^ ^- (114) 

m \ 4m 2 no / 

for plus in front of the square root. Tracking dependence on T we get that formula (I113p 

corresponds to the polarizationless solution (15B]) . 

Dipole molecules have been used for reaching of EPBEC has large electric dipole moment 

Making estimation we find that in this, most interesting case, terms proportional to 



Po play leading role. Or more precisely, we have that AttPq/uo ^> Tno and AttPq/uq 
h 2 k 2 /2m. As consequence we can simplify formula fll 12j) and find 



m \ 4m a 
for minus in front of the square root in formula (I112p and 

2 1 fh 2 k A a + A^ , 2 4ttP 2 £; 2 \ , . 

cu 2 = - — Tn A; 2 + °— . 116 

m \ 4m cr no / 

for plus in front of the square root. 

Dipole-dipole interaction is the anisotropic interaction, so we can expect anisotropy of the 
collective excitation spectrum, especially including the fact that considered equilibrium con- 
dition is also anisotropic - equilibrium electric field gives preferential direction. Nevertheless, 
dispersion dependence (11121) is isotropic. It happens due to consideration of the longitudinal 
waves, since we have used electrostatic Maxwell equations (1 1 3 [) and (11051) . Contribution of 



the transverse electric field in the wave propagation was considered in Ref. |39|, where full 
set of the Maxwell equation along with the described in this section set of QHD equation. 
It was shown that dispersion dependence becomes anisotropic, and presented here formula 
f II 1 2 ft transforms very little with replacement Po on P) cos 9, where 9 is the angle between 
direction of equilibrium electric field and direction of wave propagation. 



XI. CONCLUSION 



We have described main results obtained at development of the QHD method for quantum 
gases. Main point of the QHD method is equation derivation from many-particle Schrodinger 
equation, which explicitly models microscopic quantum particle motion. Modeling particle 
system we include the main properties of particles and corresponding inter-particle inter- 
actions. We have obtained QHD equations for the BEC with two- and three-particle short 



35 



range interaction. We derived the quantum stress tensor, which detailed analysis allows us 
to construct QHD equations. For two-particle interaction we have considered the quantum 
stress tensor up to the third order by the interaction radius, and for the TPI we have lim- 
ited our consideration by the FOIR. We have specially considered long-range interaction 
between electrical dipole moments of neutral particles being in the BEC state. In this case 
we have had to derive equations for polarization evolution in addition to the continuity and 
the Euler equation, which make up usual QHD. By means of the QHD method we consider 
influence of non-zero temperature on evolution of bosons deriving two-fluid model and con- 
sidering energy balance equation for non-condensed bosons. We have obtained contribution 
of described interactions dispersion of collective waves in the BEC. 
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